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Abstract 

Magnetocaloric effect in a double-tetrahedral chain, in which nodal lattice sites occupied by the localized Ising spins regularly 
alternate with three equivalent lattice sites available for mobile electrons, is exactly investigated by considering the one-third 
electron filling and the ferromagnetic Ising exchange interaction between the mobile electrons and their nearest Ising neighbours. 
The entropy and the magnetic Griineisen parameter, which closely relate to the magnetocaloric effect, are exactly calculated in 
order to investigate the relation between the ground-state degeneracy and the cooling efficiency of the hybrid spin-electron system 
during the adiabatic demagnetization. 

Keywords: spin-electron double-tetrahedral chain, magnetocaloric effect, entropy, magnetic Griineisen parameter, exact results 
PACS: 05.50.-Hq, 75.10.Pq, 75.30.Sg, 75.30.Kz 


1. Introduction 


The magnetocaloric effect (MCE), which is characterized by 
an adiabatic change of the temperature (or by an isothermal 
change of the entropy) under the variation of the applied mag¬ 
netic field, has a long history in cooling applications at various 
temperature regimes 111]. Since the first successful experiment 
of the adiabatic demagnetization performed in 1933 0, the 
MCE is a standard technique for achieving the extremely low 
temperatures 0. 

The theoretical prediction and description of materials with 
an enhanced or even giant MCE create real opportunities for 
the effective selection of the construction for working magnetic 
refrigeration devices. Of particular interest is the investigation 
of the MCE in various one-dimensional (ID) quantum spin sys¬ 
tems HQ or hybrid spin-electron models lldl - ll^ . The rea¬ 
son is a possibility of obtaining the exact analytical or numeri¬ 
cal results as well as a potential use of these models for the ex¬ 
planation of MCE data measured for real magnetic compounds. 
In particular, ID models may give correct quantitative descrip¬ 
tion of real three-dimensional (3D) magnets, when appropriate 
rescaling of material parameters are taken into account 111714201] . 

In the present Letter, we will investigate the MCE in a hy¬ 
brid double-tetrahedral chain composed of the localized Ising 
spins and mobile electrons, which is exactly solvable by com¬ 
bining the gen eralized decoration-iteration mapping transfor¬ 
mation 1 2j. - 241 and the transfer-matrix technique 1 3 3) ■ As 
has been shown in our previous works il5lll6|] . the considered 
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ID spin-electron system provides a suitable prototype model 
for theoretical investigation of the relation between the ground- 
state degeneracy and the cooling efficiency of the system during 
the adiabatic (de)magnetization. 

The Letter is organized as follows. In Section |2l we will 
describe the investigated spin-electron double-tetrahedral chain 
and then, the basic steps of an exact analytical treatment of the 
model will briefly be recalled. Exact calculations of the basic 
thermodynamic quantities, such as the Gibbs free energy, the 
total magnetization, the entropy and the magnetic Griineisen 
parameter, will be realized in this section. In Section [3] we 
will particularly discuss the numerical results acquired for the 
ground state, the entropy and the magnetic Griineisen parame¬ 
ter as functions of the applied magnetic field under the assump¬ 
tion of the one-third electron filling of each triangular cluster 
and the ferromagnetic exchange interaction between the mobile 
electrons and their nearest Ising neighbours. Einally, the Letter 
ends up with a summary of our findings in SectionH) 


2. Model and its exact solution 


Let us consider a magnetic system on a double-tetrahedral 
chain, where nodal lattice sites occupied by the localized Ising 
spins regularly alternate with three equivalent lattice sites avail¬ 
able to two mobile electrons. The magnetic structure of the con¬ 
sidered ID spin-electron model is schematically illustrated in 
Fig.IU Assuming the exchange interactions between the near¬ 
est neighbours, the on-site Coulomb repulsion U > 0 between 
two electrons of opposite spins occupying the same lattice site 
and the action of the external magnetic field B on the mobile 
electrons and the localized Ising spins, the total Hamiltonian of 
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Figure I: A part of the spin-electron system on the double-tetrahedral 
chain. Full circles denote nodal lattice sites occupied by the local¬ 
ized Ising spins, while the empty circles forming triangular clusters 
are available to mobile electrons. 


where /3 - l/T is the inverse temperature (we set the Boltz¬ 
mann’s constant A:b = 1) and N is the total number of the nodal 
lattice sites (the localized Ising spins). The explicit expessions 
of the mapping parameters A, J^f j and Hefj emerging in Eq. (I2]l 
can be obtained from the ’self-consistency’ condition of the ap¬ 
plied decoration-iteration transformation (see Eqs. (9) and (10) 
in Ref. ifl^ l. At this stage, the exact calculation of the partition 
function Z, of the spin-electron tetrahedral chain is formally 
completed, because the partition function Zi of the spin-1/2 
Ising linear chain in a magnetic field is known | 


the model reads; 

iUi) relT.il 




(jM) 


H, 




(1) 


In above, the summation (/, j) runs over the lattice sites form¬ 
ing triangular clusters, while the summation {j, k) runs over the 
lattice sites of triangular clusters and the nearest-neighbouring 
nodal lattice sites. The operator cj^ (c,\s) represents usual 
fermionic creation (annihilation) operator for mobile electrons 
occupying the /th lattice site with the spin s e (t, i), nj^s - 
c^j fj,s is the number operator of the mobile electron at the /th 
lattice site and cr^ labels the localized Ising spin at the kth nodal 
lattice site. The hopping parameter f > 0 takes into account the 
kinetic energy of mobile electrons delocalized over triangular 
clusters and J stands for the Ising-type coupling between mo¬ 
bile electrons and their nearest Ising neighbours. Einally, the 
last two terms in Eq. ([TJ represent the Zeeman’s energies of the 
mobile electrons and the localized Ising spins, which may be 
in general different due to a difference in the relevant Lande g- 
factors absorbed into the definition of the ’effective’ magnetic 
fields He - gefisB and Hi - gifisB (ge is Lande g-factor of 
the mobile electrons, gi is Lande g-factor of the localized Ising 
spins and fiB is Bohr magneton). 

2.1. Partition function 

One can note that the spin-electron model under consider¬ 
ation can alternatively be viewed as the spin-1/2 Ising linear 
chain, whose bonds are decorated by triangular clusters avail¬ 
able for two mobile electrons. Erom this point of view, the par¬ 
tition function of the system can exactly be derived within the 
generalized decoration-iteration mapping transformation iB- 
2411 (as described in our previous work ilhll l. As a result, one 


obtains a simple relation between the partition function Z of 
the investigated spin-electron tetrahedral chain and the partition 
function Zi of the uniform spin-1/2 Ising linear chain with the 
effective nearest-neighbour coupling Jeff and the effective mag¬ 
netic field Heff ; 

Z iP, J, t, U, Hi, He) = A^ZliP, Jeff, Heff), 


2.2. Magnetization, entropy and magnetic Griineisen parame¬ 
ter 

In this part, we present the exact solution for the Gibbs 
free energy the total magnetization M, the entropy S and 
the magnetic Griineisen prameter E// of the investigated spin- 
electron tetrahedral chain. The first three physical quantities 
immediately follow from the relation (|2]); 
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The partial derivatives of the functions InZ/ and In A, appear¬ 
ing in Eqs. (|4li and (|5]l, satisfy the general equations: 
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where i = sin\\{pHeffl2^, c - cosh(y6//ey//2) and Q — 

^sW{pHeff/2) ■+ exp(^-pjeffj. Eorms of the functions W:^, 
Wo and W are listed in Eq. (10) of Ref. ifl^ . 

The so-called magnetic Griineisen parameter T h, which can 
be calculated from the relation (see Ref. ll27n for a recent re¬ 
view): 
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(Ch is the specific heat at the constant magnetic field H), refers 
to a thermal response of magnetic system with respect to a vari¬ 
ation of the external magnetic field. This physical quantity rep¬ 
resents a ma gne tic analog of the classical thermal Griineisen pa¬ 
rameter 12814301] . It is noteworthy that the magnetic Griineisen 
parameter ([Sll diverges at a quantum phase transition driven 
by the external magnetic field quite similarly as the thermal 
Griineisen parameter does at a quantum phase transition driven 
by the external pressure. From this point of view, the magnetic 
Griineisen parameter ([8]) provides a valuable tool for an experi¬ 
mental identification of the field-tuned quantum phase transi¬ 
tions 113 114331] . In addition, it is clear from Eq. ® that the 
magnetic Griineisen parameter Fh is proportional to the adi¬ 
abatic cooling rate {dTjdH)^ and hence, it represents a key 
physical quantity for an investigation of the cooling efficiency 
during the adiabatic demagnetization especially in a vicinity of 
field-induced phase transitions. A direct substitution of the to¬ 
tal magnetization (|4]i and the temperature derivative of the en¬ 
tropy into the expression ([8]l yields the following form of 
the magnetic Griineisen parameter for the hybrid spin-electron 
double-tetrahedral chain: 
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(5) 


The second partial derivatives of the functions InZi, In A that 
emerge in Eq. (|3i can be obtained by differentiating Eqs. (|6l) 
and (|7]) with respect to the relevant variable. It should be noted, 
however, that the resulting expressions for these derivatives are 
too cumbersume to write them here explicitly. 


3. Results and discussion 

In this section, we present the most interesting results ob¬ 
tained for the spin-electron double-tetrahedral chain by consid¬ 
ering the particular case with the ferromagnetic Ising interac¬ 
tion J < 0 between the localized Ising spins and mobile elec¬ 
trons. To reduce the number of free interaction parameters, we 
will also assume equal ’effective’ magnetic fields H = Hj - Hg 
acting on magnetic particles. 


3.1. Ground state 

Eirst, let us comment on possible phases that may appear in 
the ground state of the investigated spin-electron model. A typ¬ 
ical ground-state phase diagram, constructed in the f/|7| - HI\J\ 
plane by assuming various values of the Coulomb term t//|7|, is 
displayed in Eig.|2l As one can see from this figure, the ground 
state of the system contains two possible phases, namely, the 
ferromagnetic (EM) phase and the frustrated (ERU) phase. The 
boundary, which represents the first-order transition between 
the relevant phases, is given by the condition: 


1^1 
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( 10 ) 


It is worthwhile to remark that the EM and ERU phases also 



Figure 2: The ground-state phase diagram of the spin-electron double- 
tetrahedral chain with the ferromagnetic coupling J <Q constructed in 
the tl\J\ - HI\J\ plane for the Coulomb terms Ul\J\ = 0,5,10,30. 


appear in the ground state of the antiferromagnetic counterpart 
of the model (see our preceding paper ITi^] ). Thus, we present 
here just their eigenvectors and brief definitions for the sake of 
easy reference: 

• The ferromagnetic (EM) phase: 


|EM> = f]|T)o 


k=\ 
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( 11 ) 




where oj - (i^ = -l). In this phase, the mobile 

electrons at each triangular cluster underlie a quantum 
superposition of three ferromagnetic states .|.|0), 

c^ 2|C)[3||0)> |I0) ^ttd the Ising spins localized at 

nodal lattice sites occupy the spin state cr^ = 1 /2. 

The frustrated (ERU) phase: 


|ERU> = 


niT(i)>a 


[ sin ip (c^i 


,t 


k=i 


V 6 


‘'L.i'-L.T ‘'l3,i‘'ll,T^ 


V 2 


cos^J]4.tC^J|0> (// = 0) 

j=i 


k=l 


n I [ sin (p 


'‘^/t2,i%,T' 


k3,l 
/t3,T'/tl,T 


cL .cl, J 


■ V2cos(^^ VIJIO) (^>0), (12) 

2=1 


where tanp = + 2t + sj(U + Ifp- + 32f -). In 

this phase, the mobile electrons from each trian¬ 
gular cluster show the quantum entanglement of 
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six intrinsic antiferromagnetic states 
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10 ), 

and three non-magnetic ionic states 
(j - 1,2,3), while the arrangement of the localized 
Ising spins depends on a presence of the external magnetic 
field; if // = 0, the Ising spins are completely free to flip 
in arbitrary direction, while they are fully polarized into 
the field direction, if H 


In both Eqs. (fTTI) and ([12]), the products run over all primitive 
unit cells, the state vector | T)o-t (I i)o-t) determines the up (down) 
state of the Ising spin localized at ki\\ lattice site and |0) labels 
the vacuum state. As has been evidenced in Ref. ifT^l . the above 
ground states are macroscopically degenerate, which conse¬ 
quently leads to a residual entropy 5/3A = ln203 ss 0.231 
in both the phases. The FM phase is macroscopically degen¬ 
erate due to chiral degrees of freedom of the mobile electrons, 
while the FRU phase exhibits a macroscopic degeneracy ow¬ 
ing to a kinetically-driven spin frustration of the localized Ising 
spins caused by the antiferromagnetic alignment of the mobile 
electrons. Arbitrary but non-zero magnetic field tends to align 
the Ising spins into the field direction and thus, it cancels the 
macroscopic degeneracy of the FRU phase (and also the asso¬ 
ciated residual entropy SI3N = In 2'^^). By contrast, the FM 
phase remains macroscopically degenererate in the whole pa¬ 
rameter region. 


3.2. Adiabatic (de)magnetization process 

Now, let us turn to the discussion of the magnetocaloric prop¬ 
erties of the investigated model. In Fig. |3]we depict isother¬ 
mal changes of the entropy per one magnetic particle S/3N 
(recall that the system is composed of N Ising spins and 2N 
mobile electrons) under the field variation assuming the fixed 
Coulomb term U/\J\ - 5.0, the hopping parameter t/\J\ = 1.5 
and various temperatures. As one can see from this figure, 
the entropy isotherms monotonously decrease from its maxi¬ 
mum at H/\J\ - 0 upon the increasing magnetic field down to 
temperature T /\J\ — 1.0. Below T /\J\ - 1.0, the entropy ex¬ 
hibits non-monotonous dependencies as a function of the ex¬ 
ternal magnetic field with a pronounced peak at the critical 
field ([Tol l, at which the system undergoes a phase transition be¬ 
tween the FRU and FM phases. Finally, the continuous entropy 
isotherms split into the isolated points with the coordinates 
[H/\J\,S/3N] = [0,ln2'/3], [H/\J\,S/3N] = [HJ\J\M3^/^] 
and the lines S/3N - 0 (H < He), S/3N - In 2*^^ (H > He) 
when the temperature reaches the zero value. 

In general, the residual entropy found in the frustrated mag¬ 
netic systems and at critical points corresponding to the phase 
transitions between different magnetic structures gives a rise to 
an enhanced MCF accompanied by a relatively fast cooling of 
systems during the adiabatic (de)magnetization. To investigate 
the efficiency of the adiabatic cooling of the considered spin- 
electron system in these regions, we plot in Fig. |4|the magnetic 
Griineisen parameter multiplied by the temperature TFh versus 
the external magnetic field for the relatively low temperature 
T l\J\ - 0.2 and a few values of the parameters f/|7| and t//|7|. 



Figure 3: The isothermal dependencies of the entropy 5/3A on the 
magnetic field HI\J\ for the fixed Coulomb term Ul\J\ = 5.0, the 
hopping parameter tl\J\ = 1.5 and the various temperatures Tl\J\ = 
0,0.1,0.2,..., 1.5 (from bottom to top). 


Evidently, the low-temperature TVh curves depicted in this fig¬ 
ure exhibit two local maxima at very low (but non-zero) mag¬ 
netic field and slightly above the critical field (ITOl l of the zero- 
temperature phase transition FRU-FM in addition to the one lo¬ 
cal minimum slightly below (ITOl i. The sign change of the TTh 
product observed close to the critical field (ITOl i clearly indicates 
a rapid accumulation of the entropy due to a mutual competi¬ 
tion between the neighbouring ground-state configurations (see 
e.g. the curves plotted for f/|y| = 1.5 inFig.|4^andfor f//|y| - 5 
in Fig.llJ) and compare them with the corresponding isothermal 
dependence of SI3N in Fig. [3]). Furthermore, the peaks ap¬ 
pearing around the field-induced phase transition FRU-FM are 
much higher than those observed at relatively small magnetic 
fields. For example, for t/\J\ = 1.5 and U/\J\ - 3,5,10 we 
have TTh ~ 41.6,19.8,6.5 at///|7| 2.7,2.3,1.8 and TTh ~ 

15.3,7.1,2.1 at///|7| 0.007,0.014,0.043 (see Fig. Hh). It 

is thus clear that the cooling capability of the system during 
the adiabatic (de)magnetization is substantially higher (approx¬ 
imately three times) just above the critical field, where strong 
thermal excitations of the mobile electrons are present at rel¬ 
atively low temperatures due to breaking up the quantum su¬ 
perpositions of their up-down states, than that one at relatively 
small magnetic fields, where the Zeeman’s splitting of energy 
levels of the frustrated Ising spins takes place. Except to this 
behaviour. Fig. |4| also illustrates the effect of the hopping pa¬ 
rameter and the on-site Coulomb repulsion on the enhancement 
of the MCE in the investigated model. It is quite evident from 
this figure that the adiabatic cooling rate of the system generally 
increases in the vicinity of the zero field and nearby the field- 
induced phase transition FRU-FM with increasing the kinetic 
term t/\J\ (see Fig.|4^), while it decreases in these regions with 
increasing the Coulomb parameter U/\J\ (see Fig.|33). 

To discuss the MCE in the considered spin-electron double- 
tetrahedral chain, we may alternatively investigate an adiabatic 
change of temperature of the model under the magnetic field 
variation. For this purpose, we present a density plot of the en¬ 
tropy as a function of the magnetic field and temperature for 
the fixed Coulomb term U/\J\ = 5.0 and the hopping param- 
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Figure 4: The magnetic Gruneisen parameter multiplied by the temper¬ 
ature TTu versus the magnetic field for the fixed temperature Tl\J\ = 
0.2 and the set of parameters: (a) UI\J\ = 5.0, tl\J\ = 1.0,1.3, 1.5, (b) 
tl\J\ = 1.5, UI\J\ = 3,5, 10. 


eter f/|y| = 1.5 in Fig. |5] Isentropic changes of temperature 
upon varying the magnetic field can be identified in this fig¬ 
ure as contours of the constant entropy displayed by solid and 
broken lines. In accordance with the previous discussion, the 
enhanced MCE during the adiabatic (de)magnetization can be 
found just above the zero field as well as nearby the critical field 
corresponding to the field-induced phase transition between the 
FRU and FM phases. The most abrupt drop in temperature up to 
the zero value is achieved under the adiabatic condition in these 
regions if the entropy of the system is set sufficiently close to 
the values SIW = \n2^l^ ^ 0.231 and 5/3A^ = In 3*^^ 0.366, 

respectively. 

4. Conclusions 

The present Letter deals with the ground-state and magne¬ 
tocaloric properties of the double-tetrahedral chain, where the 
nodal lattice sites occupied by the localized Ising spins regu¬ 
larly alternate with three equivalent lattice sites available for 
two mobile electrons. Based on the exact solution of the model 
presented in Ref. ns, we have analytically derived exact re¬ 
sults for the basic thermodynamic quantities, such as the Gibbs 
free energy, the total magnetization, the entropy and the mag¬ 
netic Gruneisen parameter of the system. By considering the 
ferromagnetic exchange interaction (7 < 0) between the mo¬ 
bile electrons and their nearest Ising neighbours, we have found 
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Figure 5: A density plot of the entropy as a function of the magnetic 
field and temperature by assuming the fixed Coulomb term U l\J\ = 5.0 
and the hopping term t/|7| = 1.5. The displayed curves correspond to 
isentropy lines, namely, S/3N = 0.01,0.1,0.2 (black broken curve), 
S/3N = ln2‘/3 (white broken curve), S/3N = 0.24,0.26,..., 0.4 
(black solid curves) and S/3N = In 3*^^ (white solid curve). 


two macroscopically degenerate ground states (FRU and FM), 
which can also be found in the antiferromagnetic counterpart of 
the model M- We have shown that the macroscopic degener¬ 
acy of the FRU phase arising due to the kinetically-driven frus¬ 
tration of the localized Ising spins as well as the macroscopic 
degeneracy of the system at the field-induced phase transition 
between the FRU and FM phases perfectly manifest themselves 
in the enhanced MCE during the adiabatic (de)magnetization. 
It has been evidenced that the cooling capability of the sys¬ 
tem during the adiabatic (de)magnetization is approximately 
three times higher nearby the phase transition FRU-FM, where 
strong thermal excitations of the mobile electrons are present at 
relatively low temperatures due to breaking up the quantum su¬ 
perpositions of their up-down states, than that one at relatively 
small magnetic fields, where the Zeeman’s splitting of energy 
levels of the frustrated Ising spins takes place. 
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